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^ ■ Abstract 



^ ■ We discuss observational constraints coming from CMB and type la supernovae, 

CN . for the model of accelerated universe produced by gravitational leakage into extra 

dimensions. Our fits indicate that the model is currently in agreement with the data. 
'-^1 We also give the equations governing the evolution of cosmological perturbations. 

" " Future observations will be able to severely constrain the model. 



1 Introduction 



Supernovae observations have recently provided evidence that the expansion of the 
5^ \ Universe is undergoing a late time acceleration ^. This acceleration can be 

explained in the framework of standard cosmology by a non vanishing cosmological 
constant. Although in agreement with current observations, such an explanation 
exacerbates the usual cosmological constant problem because it requires an expla- 
nation for its very small, but non zero, value. 

One may wish to find alternative explanations for the acceleration, and there are 
several proposals in the literature. Here we explore a scenario proposed in |§, ^, 
based on the model of Dvali-Gabadadze-Porrati of brane- induced gravity . This 
proposal explains the observed late time acceleration of the expansion of the Uni- 
verse through a large scale modification of gravity coming from "leakage" of gravity 
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at large scale into an extra-dimension, and without requiring a non vanishing cosmo- 
logical constant. The interesting point about this model from a phenomenological 
perspective is that it is a testable alternative to a cosmological constant model with 
the same number of parameters. This is in contrast to models of "quintessence" 
where the equation of state of the new component becomes a free function that 
needs to be constrained. 

In 1^ it has been shown that the model was in qualitative agreement with all 
known cosmological observations. The purpose of this work is to go one step fur- 
ther and quantitatively confront the model with observations of supernovae and the 
cosmic microwave background (CMB). 

The paper is organized as follows, in section |^ we discuss the dynamics of the 
background metric of the universe in the model. We first introduce in a few words 
the brane- induced gravity model of Dvali-Gabadadze-Porrati (see also [0, P]) 
which provides the framework (subsection |2.1| ). We then discuss the cosmological 
dynamics for the accelerated solution considered in this paper (subsection p.2|) . In 
the following, we confront the model with the Supernovae observations of the Su- 
pernova Cosmology Project (SCP) M (subsection lO) and CMB data (subsection 



3J). 

Our fits indicate that the model is currently in agreement with SNIa and small 
scale CMB data. One can hope to discriminate the model from standard cosmol- 
ogy using future precision cosmological parameters measurements, but also maybe 
modifications in the growth of large scale structures. 



2 Model definition and Background dynamics 

In the following subsections we summarize the main features of the model under 
consideration, study the dynamics of the background metric and confront the pre- 
dictions of the model to the supernovae observations. 

2.1 Brane-Induced Gravity Models in a Few Words 

The brane-induced gravity models are a particular class of brane-world models, 
which can be defined as models where our four dimensional (4D) universe is con- 
sidered to be a surface (called brane) embedded into a higher dimensional bulk 
space-time. 

Brane-world models are inspired by superstring-M theory, and can be regarded 
as some low energy effective models of more fundamental underlying theories, but 
have also interest on their own in providing new phenomenological ideas. We will 
only consider here the case where the bulk is five dimensional (5D). The brane em- 
bedding into the bulk is defined by the coordinates X^{x^) of the brane worldvolume 
(parametrized by coordinates x^) into the 5D space-time. The dynamics of gravity 
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is governed by the usual 5D Einstein-Hilbert action 



Seh = ^ f d'X J\^\ (5)/?, (1) 



where M(5) denotes the 5D reduced Planck mass. The bulk metric ^^^qab induces 
through the embedding X^{x^) a metric g^j^^ on the brane (called induced metric) 
defined by^. 

g^u =('^ gABd.X^'d.X^. (2) 

In the above equation, we have put an upper index "^^^ on quantities (e.g. the 5D Ricci 
scalar '^^^R or the 5D metric ^^''gAB) to distinguish them from their 4D counterparts 
depending only on the induced metric (e.g. R or (yf^jy). 

In the brane- induced gravity models 0, 0, ||, ^ [1^, the gravitational action 
contains an other term Seh-, in addition to the 5D Einstein-Hilbert term (||), given 
by 



Seh = ^ dS ^\g\ R. (3) 

Z J brane 

This term is the usual 4D Einstein-Hilbert term computed here on the brane and 
with the induced metric, with Mpi a mass parameter. The latter can be interpreted 
as the usual 4D reduced Planck mass, from the calculation (see below) of the force 
between two static massive sources on a flat brane and bulk background!]. The 
origin of Seh in brane world models is discussed in more details in [|, 0, §|. It arises 
generically from quantum correction coming from the coupling between bulk space- 
time and brane localized matter flelds when the conformal invariance of the brane 
theory is broken (see e.g. fl^). In the model at hand the dynamics of gravity is then 
governed by the sum of the two kinetic terms Seh and Seh- 

As a consequence of the presence of the brane-induced term (^), one can show 
than the gravitational force experienced by two static point like sources located 
on the brane is the usual 4D gravitational force for distances smaller than the 
crossover scale deflned by 

- - (4) 



For distances larger than Tc, on the other hand, the force turns to a 5D regime 
where it follows the 5D behavior. On scales smaller than M'^^ one also expects 

^In the following, we use upper case Latin letters A,B, ... to denote 5D indices, Greek letters 
fj,,^,... to denote indices parallel to the brane world volume, 5 an index transverse to the brane, 
and Latin letters ... to denote space-like indices parallel to the brane world volume 

^We will not address here the issue of the vDVZ discontinuity, see and jl2 , 13 for discussions 
of this issue. 



3 



modifications in the gravity law, however for the parameter choice relevant to this 
work, the modifications occur on scales much smaller than those accessible by gravity 
experiments 0. 

This perturbative behavior has an exact parallel in cosmology, where one can 
show 1^ that, for a Z2 symmetric brane world (see []T3|, [T^ ^ for discussions of 
cases where the Z2 symmetry is relaxed), the expansion of the Universe is governed 
by the usual 4D Friedmann's equations whenever the Hubble radius is smaller 
than Tc, and enters into a non conventional regime for larger Hubble radii. 

In the following subsection we will discuss in greater detail this cosmological 
evolution. At this point let us first say that an obvious criteria that the model should 
reach in order to comply with the known behavior of gravity at large distance, as well 
as with the observed cosmology, is that should be made large enough. The more 
stringent limit comes indeed from cosmology requiring Tc to be of the order of, or 
larger than, the today's Hubble radius . When Tc ~ ■, one thus expects that 
cosmology is very close to standard cosmology up to very late time, and in particular 
all successes of standard cosmology such as BBN are left unchanged by this choice of 
parameters. However the very recent evolution of the universe is different. Indeed, 
as will be reminded in more detail in the next subsection, a particular class of 
solutions shows a late time accelerated expansion without the need for a non zero 
cosmological constant. For values of Tc of order as needed to fit the Supernovae 
observations (see section |3.1| ), one finds from Eq. (^) that M(5) ~ 10 — 100 MeV. 
Such a low value of the 5D Planck mass is perfectly consistent with observations 
and high energy experiments as shown in |^. Induced-gravity models have been 
shown to provide a framework for realizing models with a very low quantum gravity 
scale without confiicting with any experimental facts P]. 

2.2 Background cosmological dynamics 

In the model considered here, the geometry of our 4D Universe is at all time described 
by an ordinary FLRW space-time with a line element of the form 



where -0 are angular coordinates, k = —1, 0, 1 parametrizes the brane world spatial 
curvature, and Sk is given by 



The cosmological standard observers are assumed, as usual, to be at rest with respect 
to the comoving coordinates x*. The only difference with standard cosmology is in 



Q jj^j^dx doc 
— dt^ + a^{t) dx^ dx-' jij 



(5) 
(6) 
(7) 
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the dynamics of the metric which is encoded into Friedmann-hke equations different 
from to the ordinary 4D ones. For a given content of the universe, with total energy 
density p (and pressure p) , the standard first Friedmann's equation is now replaced 
by i 
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where 



H.i% (10) 

a at 

is the Hubble parameter of our universe^. The energy-momentum conservation 
equation, on the other hand, takes the usual form 

p + 3H{p + p) = 0. (11) 

Equations (^ and (|TT]) are all what is needed to characterize the cosmology we are 
interested in here. They lead to. 



H\z) = Hi { 



^^(1 + ^)2+ a^ + ^fi,(l + z)3(i+-.) , (12) 



where z is the redshift and we have assumed that p is given by the sum of the energy 
densities pa of different components (labeled by a) with constant equation of state 
parameters Wa- The fis for matter and curvature are defined in the usual way by 



-k 



^ (14) 

(15) 



whereas flr^ is given by 



The normalization condition for the Qs, 



+ J^r. + , ^r. +J2nJ =1, (17) 



* There is another set of solutions for a Z2 symmetrical brane. Those were derived in ^ and 
are not considered here. 
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differs from the usual relation fi^ + J2a = 1- 

Equation (j^) implies that whenever p/Mp^ is large compared to (or in other 
words, whenever is small with respect to r^), the cosmological evolution fol- 
lows that of standard cosmology. In this case equation (|^) reduces to the standard 
Friedmann's equation 



2 , k p 



a2 SM^i' 



When (and if) p is driven to smaller values by the cosmic expansion, the expansion 
of the Universe enters into a non conventional phase and asymptotes to a de Sitter 
solution when p becomes negligible w.r.t. Mp^/r^. One has a transition to an 
accelerated expansion happening approximately when the Hubble radius crosses 
the threshold Vc- We would like to stress that this last accelerated phase is not 
triggered by a cosmological constant (that can be consistently set to zero) but is 
due to the presence of two kinetic terms for the graviton in the action. Namely, 
bulk gravity sees the induced kinetic term on the brane (^ as a source term, and 
for an empty universe, there is a self-inflationary solution]^ to Einstein's equations 
to which a universe with decreasing energy will asymptote. This solution acts as a 
late time attractor to early standard cosmology. 

In the following we will then only consider a universe with a zero cosmological 
constant, and usual (dark, baryonic, ...) matter content. One can further notice that 
the above described cosmology is also exactly reproduced by standard cosmology 
with a dark energy component with a 2;-dependent equation of state parameter 
Wx^{z). For a universe containing only non relativistic matter, (z) is given by 
(see i) 

= .... w / 1 I , - 1- (19) 



At large redshift tends toward —1/2 reflecting the fact that the dominant term 
in equation (0), after matter and curvature terms, redshift as (1 + z)^^"^ at large z. 
At low z, however, Wx^ decreases toward an (0^, nM)-dependent asymptotic value. 
For a fiat universe, the latter is simply given by|^ — 1/(1 + Qm)- 

In the following sections we give the results of fitting SNIa and CMB observables 
with different cosmological parameters in the framework of the cosmology defined 
by equations (j^) and ([Tl|). We will denote a set of cosmological parameters such 
as flr^ or flM characterizing a given cosmology. 



^this solution is in a way the late time analog of Starobinsky's first model of inflation where 
terms quadratic in the Ricci tensor are sourcing similarly a self-inflationary solution [ p^ . 
^°e.g. for = 0.3 and fc = 0, w'^^ at low z tends toward —0.77. 
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3 Confrontation with observations 



3.1 Confrontation with supernovae observations 

We have fitted the supernovae data set from the SCP with the luminosity distance 
calculated using Equation (|l^). Because the geometry of the Universe is given by 
usual FLRW (|) one can use the standard formula for the luminosity distance cIl as 
a function of the redshift z, 



.,^^,-■(1 + ^) ^^^^"''') . (20) 



with dc{z) defined by 



z 



and H{z) given by equation (|12|). We then use this luminosity distance to fit the 
data. The fit is done using 4 free parameters: the cosmological parameters = 
{QMy^rJ, the intrinsic magnitude of the supernova A4, and a parameter a related 
to the intrinsic luminosity-decline rate relation (stretch factor s). is given by 

^.2(g ^ _ ^ (7W + a(l - s,) + 5 logioidije, Zj)) - mQ^ ^^2) 

i=l 

The data set consisting of 54 supernovae (18 nearby ones and 36 at high redshift) 
is shown in figure 0. Since we assume no prior knowledge of the parameters and as 
we are not interested by a and A^, we have to marginalize over them. We do this in 
a Bayesian framework assuming fiat priors and Gaussian errors. These integrations 
can be carried out analytically, as shown in [^. We quote the results in appendix 
^ We have then computed confidence contours for the models 6 = {Qm, Qr-J model 
with no prior on the cosmology. These contours are plotted in fig ^ 

Assuming a spatially flat space-time, one is only left with one free parameter 
(after integration over Ai and a), e.g. Qm- is then given by the normalization 
condition (|T^ 

= ) ,nr,<l and VIm < 1- (23) 



2 

The results of the minimization gives for a fiat universe (one sigma levels) 

Qa/ = 0.181°:°^ or a, = 0.17+°:°^^ (24) 

with = 57.96, for 52 (54 SNe - 2 parameters) degrees of freedom |^ This best fit 
model is shown in figure |l|. Equation (p^ leads to an estimate for Tc in terms of the 

These numerical results are in agreement with the fit done in [ pO| . We however disagree with 
the conclusions of that work as will be discussed later (see also Ref. for a discussion of this 
paper). Note in particular that, contrarily to the claims made in [ pO[ , the z = 1.7 supernova of 
Ref. 1^, is fitted as well by the model considered in this paper, or by standard cosmology with a 
cosmological constant (see figure |^). 
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Figure 1: Magnitude vs redshift diagram for the SNIa data of reference 0] used in 
this paper. All magnitude are plotted respectively to an empty universe {Qm = 
and f^A = 0). Over-plotted are three different flat cosmological models: the best fit 
flat model in standard cosmology (with Qm = 0.28 and Q\ = 0.72, solid line), in 
the gravitational leakage cosmology (with Qm = 0.18, dotted line) and a fiat model 
in the gravitational leakage cosmology with Qm = 0.3 (dashed line). We also show 
two approximate confidence level interval for the z = 1.7 supernova of Ref. [Q], the 
outer light-gray surface represents the 95% confidence interval, the inner dark gray 
surface represents approximately the 68% confidence interval. This last supernova 
was not included in the fit. The values of a (related to the stretch factor) and Ai 
(intrinsic magnitude) have been fitted independently for all the models. The data 
are plotted here with a = 0.6. 
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Figure 2: 68.3%, 90% and 99% confidence regions for (fiAf,^rJ i^i gravitational 
leakage scenario, assuming no priors knowledge of a and M.. 

Hubble radius Hq"^ given by 

re = \.2\^^^E^\ (25) 
3.2 Confrontation with CMB Observations 

Another set of cosmological observables which has recently been measured with great 
precision is the CMB temperature power spectrum. In this subsection we would like 
to compare the predictions of the model considered in this paper to the results of 
these observations. 

For this purpose, we used a modified CMBFAST replacing the first Fried- 
mann's equation by equation (||). The equations for the growth of cosmological 
perturbations were kept the same as in usual cosmology (except for the background 
evolution). As is discussed qualitatively in appendix B, this is justified for the small 
scale perturbations and for processes happening early enough in the history of the 
Universe. On the other hand, one can expect deviations from the standard picture 
at large scale (and late time) where (and when) the effect of the extra dimension 
began to be felt. This concerns scales of order of today's Hubble radius and pro- 
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cesses happening in the late history of the universe. A more refined discussion of 
this, which involves the integration of bulk equations of motions for perturbations, 
is left for future work |^ . 



We explored the six-dimensional parameter space, 6 = (fit, ^2^^; i^d, i^b, 
where Ud = flcdmh"^, = ^hh"^ and A and n are the amplitude and slope of the 
primordial spectrum of perturbations. We used a Markov chain method to explore 
the likelihood in this parameter space. When it has converged the method produces 
a chain of models that are sampled from the probability distribution of 6. The 
details of our procedure are given in appendix |p. 

Figure |^ shows the probability distribution for each of the six parameters ob- 
tained. As expected the CMB data prefers spatially fiat models. Figure H shows 
the results of our analysis in the Qm — plane. The shaded region was drawn to 
contain approximately 95 % of the models in our chain, the line marks the location 
of spatially fiat models. The constraint on is coming mainly from the position 
of the acoustic peaks so there is a natural degeneracy in the ^Im — plane which 
is apparent in the plot. 




Figure 3: Marginal distribution for each of the 6 parameters used. 

The probability distribution for Qk shown in figure ^ peaks around = 0, a 
spatially fiat universe. Thus it is natural to further restrict ourselves to fiat universes 
which we can do by considering only samples in our chain with negligible curvature. 
The probability distribution for Qm under this assumption is shown in figure ^ 

Figure |^ shows that a model with Qm = 0.3 provides a good fit to both SN and 
CMB data. It should be noted however that the CMB prefers a slightly larger value 
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Figure 4: Allowed region in the Qm^^vc plane (shaded). The line shows the location 
of spatially flat models. The shaded region was drawn to contain approximately 95 
% of the models in our chain. 



of than the SN, although both ranges overlap at la. In turn, the CMB can 
constrain the physical densities in matter quite accurately, uu + ~ 0.12. This 
constraint translates Qm = 0.3 to a Hubble constant h ^ 0.63 in good agreement 
with direct measurements, e.g. h = 0.72 ± 0.08 from the HST key project |^ . 



In figurey we show what we could call our "concordance" model, 6 = (^2^, ilrc^^d, 
cjb, n, A) = (0, 0.1225, 0.1, 0.02, 0.96, 0.57) which has = 0.3 and ~ 140 for the 
full data set (135 data points). For reference we also show the power spectra for a 
standard cosmology model with fl\ = 0.7, a model that has all other parameters the 
same but where the acceleration is produced by a cosmological constant. Both mod- 
els provide an acceptable fit to the current data, but because there angular diameter 
distance to the last scattering surface differs by approximately 4% they should be 
easily distinguishable by future generation of CMB experiments. The difference at 
low multipoles should be regarded with care because on this very large scales the 
physical effects we ignored could be relevant (see appendix 0) . 



4 Conclusions 

The fits done in this work show that the model of accelerated universe through 
gravitational leakage into extra dimension of Ref. ^ is in current agreement with 
SNIa and CMB data. The degeneracies in parameters estimations using one data 
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Figure 5: Marginal distributions for Q,m under the assumption that the universe is 
spatially flat. The solid line shows the results from CMB and the dashed line from 
the SN. 

set (e.g. CMB) can be partially lifted using the other (e.g. SNIa) as in standard cos- 
mology. The Supernovae data prefer a slightly lower value of (^m = 0.18l'lo;Qg) 
than the CMB for a flat universe, however a concordance model with {Qk,^rc:^d, 
uJb,n,A) = (0,0.1225,0.1,0.02,0.96,0.57) which has = 0.3 (and ~ 140 for 
the full data set (135 data points)) provide a good fit to both sets, all the more 
as we have not included systematic errors in our parameter estimations. For this 
model the crossover distance between 4D and 5D gravity is given by Tc ~ 1.4 Hq^. 

We have also given the equation of evolution for cosmological perturbations. 
Those equations were used to justify the approximation we made to compute cos- 
mological perturbations, namely we used standard four dimensional evolutions equa- 
tions over a background with a scale factor given by the accelerated solution given 
in 0. This is justified for small scale CMB anistotropies (scale smaller than the 
crossover scale Tc). From those equations, and the known behavior of gravity in 
the model at hand, one can also expect modifications in the growth of large scale 
structure. This could potentially lead to a way to discriminate between standard 
cosmology and the model considered in this work, and is left for future investigation. 

We want to end by noting that the model under consideration is very predictive 
in the sense that future observations have the potential to rule it out. In contrast 
to quintessence models, this model has the same number of free parameters as the 
usual LCDM model. With the advent of new precision cosmological measurements 
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Figure 6: Model predictions and current CMB data. The solid curve curve is for 
a model with 6 = (1]^, Wd, ujy„n,A) = (0,0.1225,0.1,0.02,0.96,0.57) while the 
dash curve is for standard cosmology with a cosmological constant I^a = 0.7 (other 
parameters were kept the same). 



such as new SNIa observations, CMB measurements, ongoing galaxy surveys such 
as Sloan and 2dF, weak lensing surveys, etc. it should be possible to test the model 
very accurately (for a recent summary of how different observations will constrain 



the matter content of the universe see Ell] and references therein) 
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A Marginalization for Supernovae 



Following |T^ the defined by equation (|2^) can be integrated analytically over 
}A and a to yield 



-2 In 



/oo / \ 

^ daexp [--XM-inti^, 



a 



(26) 
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where XM-mti^y^^) is defined by 

r / 1 

X^_int(f,a) = -2lny^^dMexpi^--x\e,M,a] 



(34) 



B Dynamics of scalar cosmological perturbations 

We briefiy summarize here the equations governing the cosmological perturbations 
in the model at hand. These equations will be derived and discussed in more details 
elsewhere, and are only given here for the case of a fiat universe. 



Our starting point is an equation derived in |2J], relating the 4D Einstein's tensor 
Gf^i, to a tensor 11^,^ quadratic in whatever source T^,^ of 5D Einstein's equations is 
localized on the brane, and a traceless tensor S^i, defined in terms of the 5D bulk 
Weyl tensor. The corresponding equation reads 



G 



1 



(35) 



with Unu given by 



n 



1 _ _ 1 _ _ 1 _ _ 1 _ 
__T T'^ -I —TT -I- — T T°'f^n —T'^n 



(36) 



and Sni, is defined by 



c _ ^5 



(37) 
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from the bulk Weyl's tensor^ bcd- '^^se T^^, is given by 

T^u = — Mp^G^u, (38) 

where T^^i, is the brane energy momentum tensor and G^i, is the 4D Einstein's tensor. 
T^i, is conserved with respect to the 4D metric on the brane, so that one has 

D^Tj: = 0, (39) 
D.fi: = 0, (40) 

where Z)^ denotes the covariant derivative compatible with the 4D metric on the 
brane, and the last equality follows from Bianchi identities. Equations (^5|) and (pOf ) 
lead to the background equation of motion ([T2|) , once one knows the background 
expression for Sfj^^. In the cosmo logical case, S^^ is in general given by some version 



of Birkhoff's theorem p5| , We have assumed for simplicity in (|T2| 



that it vanishes in the background, in which case the 5 dimensional space-time is 
simply a Minkowski space-time]^. 

We now derive from equation (|35D , the evolution equations for the cosmological 
perturbations. We write 

= ""G^^ + SG^,, (41) 

rpi. ^ ^Tl^ + STl^, (42) 

8!^ = 68!^, (43) 

where the superscript ^ stands for the background value of the corresponding tensor 
component. We define then the scalar perturbations in energy density, 6p, momen- 
tum, Sq, pressure, 6P, and anisotropic stress, Sn, for ordinary matter as 

6T^ = -6p, (44) 
57^° = V,6q, (45) 

6T; = 5P5] + {w'Wj - Sn, (46) 

where Vj is the covariant derivative adapted to the background spatial metric 
parallel to the brane. We also define similar quantities for the Weyl's fluid, following 



% 33 



= TJT^Pe. (47) 



= (48) 
58; = -^^{sPe5]+(v^V,~ls;V')sne). (49) 
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we have chosen here iniphcitly a Gaussian normal coordinate with respect to the brane 
"'^^As far as the background is concerned a non vanishing is manifesting itself as a radiation 
component into the Friedmann's equations; see Ref. where the background equations are given 
in full generality. 
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The Weyl's fluid is related to the perturbation of the bulk Weyl's tensor (gravita- 
tional waves in the bulk) through equation (p7|). 

Other useful quantities are the trace, SGt, and traceless traceless part, SGtf, of 
6G^j defined by 

SG) = 5Gt5] + (v'Vj - 5Gtf. (50) 

After some algebra, one gets then from equations (|35|) the perturbed Einstein's 
tensors over the background (pi])-([T^) 

SGtf(i ^ I = Tirf^T^-^TT^ ^ 1, (53) 

These equations replace the perturbed Einstein's equations of ordinary cosmology 
(see e.g. 



341 for a review) 



One can derive from (^) the usual conservation equations for the matter pertur- 
bations. As far as the Weyl's fluid is concerned, by taking the covariant derivative of 
equation ( |35D and using ( pQ] ) one can show that the Weyl's fluid energy density 5p£ 
is conserved but that the Weyl's fluid momentum 6q£ in general is not [^|. More- 



over one does not have an evolution equation for the Weyl's fluid anisotropic stress 
5Ti£. This means that the system of equations for cosmological perturbations do not 
close on the brane, and one needs to solve the equations of motion for gravitational 
waves in the bulk (see |^). On large scales, however, the usual adiabatic curvature 



perturbation on hypersurfaces of uniform (ordinary or Weyl) matter density is con- 
served |^T[ , since it is a mere consequence of the conservation of the energy density 
perturbation However one still cannot compute the Sachs Wolf effect because 



of the lack of knowledge of dug |^ . 

Let us make here some simple remarks. In the formalism used so far, the devia- 
tion from usual 4D cosmological perturbations equations can be separated into two 
different parts. 

We first note that the direct coupling between ordinary matter and gravitational 
perturbations is H dependent, for example one can rewrite equation (pTF) as 



.-GS = -^L-^) (55) 



with the effective direct gravitational coupling between matter and gravity given by 

1 

2H7c 



Mil = ( 1 - -i- ) . (56) 
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One can check that this couphng is never negative for the late time accelerated 
solution considered in this paper, since one always has > 1. Moreover in the 
early time of the Universe (whenever Hr^ ^ 1), one has Mp\ ~ Mp\ so that one 
can consistently ignore this effect at least up to last scattering (in contrast to what 
is happening in usual brane cosmology), which is all what matters as far as CMB 
is concerned. At the epoch of last scattering, for example, Mpi coincides with Mpi 
within a part per thousand. However, Hrc becomes of order unity at late time 
(see equation (pi])), and one can be concerned that this can have dramatic effects 
on large scale structure formation^ To be consistent one should also consider in 
this regime the effects of the Weyl's fluid source terms in the left hand side of the 
perturbed Einstein's equations ( pl| - |5^) (as well as possible non linear corrections, see 
]TT[ ) , and this can only be done properly solving for the bulk equations of motions 
for perturbations. With the formalism used so far, those source terms are the other 
manifestation of the extra dimension that we would like to discuss now qualitatively 
as far as the CMB is concerned. 

One note, that those source terms are suppressed with respect to their ordinary 
matter counterparts by a factor ifr^, in contrast to what happens usually for brane 
cosmology (see e.g @). This support the fact that when Hvc oo the theory 
looks more and more 4 dimensional. One can then start with initial conditions 
for cosmological perturbations, say after inflation, which are the one provided by 
standard 4D cosmology, and set initially all the Weyl's perturbations to zero. The 
brane perturbations will then feed up non zero perturbations in the bulk, which 
will then backreact on the brane through the Weyl's fluid perturbations leading to 
"gravitational leakage" into the extra dimension. The time scale for this leakage to 
occurs is however of order Vc which is much larger than the age of the universe at 
recombination^. 

This discussion indicates that one can consistently use the usual 4D cosmologi- 
cal perturbations equations for dealing with the growth of small scale fluctuations 
observed in CMB. The effects of gravitational leakage is then only contained in the 
background evolution, which affects the growth of the perturbations, but also the 
way they appear on the sky through a different angular diameter distance. We only 
expect possible deviations on large scale coming from the effects mentioned above, 
and also possible modifications once compared with large scale structure data. We 
let these interesting questions for future investigations, as well as a more careful 
check of the approximations made here. 



^"'This could also have potentially observable signature through standard tests of gravity. We 
however expect that when one looks at fluctuations over a given background, a local curvature 
scale should typically replace in the above equations so that we do not expect that the I 
dependence of Mp\ could have observable effects on systems where the curvature is much greater 
than today's Hq (~ r^^). This issue is likely to be related to the disappearance of the vDVZ 
discontinuity (see |l^, |l^) and will be discussed elsewhere. 

^^This qualitative picture is supported by numerical calculations in a scalar field toy model [EgI. 
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C CMB Likelihood calculation 



Here we describe the details of the CMB hkehhood calculation. 

To accelerate the calculation of the model predictions we used the k-split ap- 



proximation described in [3/]. In this approximation the high / power spectra is 
calculated in a flat model with no dark energy and then shifted appropriately in 
/ using the angular diameter distance to recombination. The likelihood for each 
model was calculated using the RADPACK]^ package. We used all currently avail- 
able CMB data (for a description of the compilation we refer the reader to the 
RADPACK documentation and to pp. 

For the present study we did not use a grid based method for calculating like- 



lihoods (such as the one described in ||3^). Following |^ we instead chose to use 



the Metropolis Hastings algorithm to generate a Markov chain of models. Because 
we are not interested in investigating multiple priors, our parameter space is rather 
small so we do not need to exploit the CMB degeneracies and we don't want to build 
a database of models to be used in future studies, the Markov-chain technique was 
very efficient and extremely easy to implement. 

In the Metropolis Hastings algorithm a chain of models is generated. Models 
are added to the chain sequentially. To find a new model for the chain values of 
the parameters are chosen at random (we choose to select models from a Gaussian 
distribution centered in the last model of the chain with a covariance matrix that is 
estimated from the chain itself). The likelihood of this new model is compared to 
the likelihood of the last model in the chain. The new model is always accepted into 
the chain if the likelihood is larger than that of the last model, if this is not the case 
it will be accepted with a probability given by the likelihood ratio of the two models. 
When the chain has converged {i. e. it has run for a sufficiently long time) each model 
can be taken as an independent sample from the probability distribution P{6\d), 
the probability of some particular value of the parameters (6) given the observed 
data (d). 

Once we have the chain of models, using histograms we can construct the distri- 
bution function of individual parameters. 
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RADPACK is publicly availa ble software package developed by Lloyd Knox. It can be obtained 



from bttp://bubba. ucdavis.edu/ 'knox/radpack.html 
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